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Abstract 

Dynamic buckling behavior of a column (rod, beam) under constant rate com¬ 
pression is considered. The buckling is caused by prescribed motion of column 
ends toward each other with constant velocity. Simple model with one degree 
of freedom simulating static and dynamic buckling of a column is derived. In 
the case of small initial disturbances the model yields simple analytical depen¬ 
dencies between the main parameters of the problem: critical force, compression 
rate, and initial disturbance. It is shown that the time required for buckling is 
inversely proportional to cubic root of compression velocity and logarithmically 
depends on the initial disturbance. Analytical expression for critical buckling 
force as a function of compression velocity is derived. It is shown that in a range 
of compression rates typical for laboratory experiments the dependence is ac¬ 
curately approximated by a power law with exponent close to 2/3. Theoretical 
findings are supported by available results of laboratory experiments. 

Keywords: dynamic buckling, Hoff problem, column, Airy equation, Euler 
force. 


1 Introduction 

Buckling of columns (rods, beams) under compression is a classical problem for me¬ 
chanics of solids. In 1744 Leonard Euler predicted critical buckling force for compressed 
column in statics. Numerous experimental and theoretical studies have revealed that 
behavior of a column in dynamics is significantly more complicated. In particular, in 
dynamics the maximum force acting in the column, further refereed to as critical force, 
exceeds Euler force. Dynamic buckling behavior significantly depends on the way of 
compression. A review of different loading conditions may be found, for example, in a 
review paper [lj. Sudden application of a force was investigated, for example, in recent 
works Mass falling on the rod was studied theoretically and experimentally in 

papers HE]. In laboratory experiments the load-bearing capacity of columns is usually 
measured in hydraulic testing machines. In this case column ends move toward each 
other with prescribed velocity mm- This loading regime is also typical for computer 
simulation of buckling at macro [8j as well as micro and nano scale levels [9]. In 1951 
Hoff has proposed the following simplified statement of this problem |B]. Compression 
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of a column with initial imperfection in a form of a sine wave was considered. It was 
assumed that the deflection of a column has the shape of the first buckling mode. This 
assumption yields the nonlinear differential equation describing buckling of a column. 
In particular, Hoff has shown that the critical force strongly depends on compression 
rate and initial imperfection of a column. Buckling under constant rate compression 
in Hoff’s statement was studied in papers The influence of axial 

inertia HU, random imperfection JT3] and boundary conditions [8] was investigated. 
Approximate solutions of Hoff’s equation are discussed in papers m HU- The de¬ 
pendence of critical buckling force on the rate of compression is obtained by fitting 
the numerical numerical results [8]. Recent experiments [7] has revealed that the de¬ 
pendence is accurately approximated by the power law. However to our knowledge 
theoretical explanation of this fact is not present in the literature. 

In the present paper simple one-dimensional model of dynamic buckling under con¬ 
stant rate compression is presented. It is assumed that the equilibrium configuration 
of the column is perfectly straight. Understanding of buckling of perfect columns is 
especially important at nano scale pun!, because nano-objects may have no defects. 
The disturbance is introduced by non-zero initial deflection. At nano scale it might 
be associated with thermal motion. Analytical dependencies of the time required for 
buckling and critical force on compression rate are derived. Theoretical findings are 
supported by the results of laboratory experiments [7]. 

2 One-dimensional model, analytical solution 

Consider simple structural model reflecting the main physical features of a column 
subjected to constant rate compression. The column is simulated by a particle (point 
mass) connected with two walls by linear springs with stiffness Cl and equilibrium 
length ao (see figure [I]). Transverse stiffness of the column is simulated by a spring 



Figure 1: Structural model for dynamic buckling of a column under constant rate 
compression. 


with stiffness Ct connecting the particle with middle point between the walls (see 
figure [I]). The walls move toward each other with constant velocity v. Potential energy 
of the system has the form: 

U = c L (Va 2 + x 2 - a 0 ) 2 + 7 } C t% 2 , (1) 

where x is a coordinate of the particle; a, a 0 are the current and initial half distances 
between the walls. 
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For any distance between the walls 2 a the straight configuration of the system x = 
0 is an equilibrium. At some distance between the walls this equilibrium becomes 
unstable and two additional stable equilibria exist: 


Xl,2 
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( 2 ) 


Here a E is a bifurcation point corresponding to critical half-distance between the walls. 
It can be show that for a > a E the straight configuration on the system (x = 0) is 
unstable. Stable equilibria are defined by formula (J2]). Corresponding Euler-like critical 
force in statics is 


Z7 i a E c T 

V E = ClCIo£e, £e = 1 -= ----, 

a 0 2 c L + c T 


( 3 ) 


where is a critical deformation in statics. Thus static behavior of the system is 
qualitatively similar to the behavior of the column under compression. Below critical 
deformation the straight configuration of the system is stable. For higher deformations 
there are two symmetric equilibria pj). Critical value of the force corresponding to 
static buckling is given by formula pL 

Consider the influence of dynamic effects on critical buckling force. Only transverse 
motions of the particle are considered. The effect of axial inertia is small and can be 
neglected HU. The equation of transverse motion of the particle is 


( 2c E CiQ \ . . . . 

mx = - 2cl + c t - , „ x, x{-t E ) = x 0 , x{-t E ) = 0, a = a E - vt. 

V V a 2 + x z J 

( 4 ) 

where m is mass of the particle, x$ is an initial deflection. Loading is applied at t — —t E: 
where t E = (ao — a E )/v is a time required for reaching Euler static force. Initial 
disturbance is simulated by non-zero initial deflection xq. Note that the statement 
used by Hoff is different [BJ- In paper [6] it is assumed that natural configuration of 
the column has finite curvature. In the present paper perfectly straight column is 
considered. 

Horizontal component of the force (see fig. [Tj) acting on the walls has the form: 


The force F has maximum value at moment of time £*, further refereed to as time 
required for buckling [ 3 ]- Parameter £* is defined by the equation: 

^F(x(t,v,x 0 ),vt)\ t=u = 0 => t*=t*(v,x 0 ). (6) 

Equation (|6| yields the dependence of i* on v in the implicit form. Substituting x(t*) 
and x 0 ) into equation (J5]) yields the expression for the critical buckling force: 


F* 


c L (a E - vt*) 


a 0 


VT a E ~ vt*) 2 + xl 



x* = x(t*, V, Xq). 


( 7 ) 
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The dependencies t*(v,Xo), x(t*,v,x o) are given in the implicit form by nonlinear 
equations Q, (|6]). 

The exact solution of equations Q, (j6]) in the general case is not straightforward. 
Therefore the following assumptions are used 

vU <C a E , |x*| <§C a E . (8) 


Linearize the expression for critical force with respect to small parameters x* and uf*: 

„2 


F* « F e + c L vU ~ c L 


aox; 
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(9) 


The relation between f* and v is derived as follows. The force F is expanded into series 
with respect to vt and x and then substituted into 


dF 


a 0 


~fa\t=t* ~ c L ( v -) = 0. 


( 10 ) 


Critical values of displacement x(£*) and velocity x(£*) are calculated using the equation 
of motion Q. The latter is linearized under assumptions 

2 c L a 0 


x 


avtx = 0, 


a = 
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( 11 ) 


The equation (11) is usually referred to as Airy’s equation (see e.g. [15]). The solution 
of equation (11) is a linear combination of Airy’s functions Ai and Bi [15] : 

x{t) = CiAi (ta^vz^J + C 2 Bi , te = a? 3 (a E — ao)u _3 , 

x 0 Bi '(te) „ xqA\'{t e ) 


C\ = 


Bi , (r i? )Ai(r £ ;) - Bi(T E )AV(T E ) ’ 


C 2 = 


Bi / (T E )Ai(r i? ) - Bi(T E )Ai'(T E )' 


( 12 ) 

For positive arguments Ai exponentially decays and Bi exponentially grows [15]. In 
particular, for large values of the argument the following asymptotic formulas can be 
used: 
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Therefore the first term in the expression (12) for x(t) is neglected. Then asymptotic 
formulas (13) and formula (10) yields the following implicit relation between £* and v: 

(14) 


vr 2 Ai r E 2 f x 0 \ 

- -rr- — Bl ( r *) Bl ( r *) > 

( 1-£ E ) 3 \ a o ) 


1 1 

r* = h« 3 ?' 3 , 


where v s = ao^/2x]~Jm is a characteristic velocity, corresponding to velocity of longi¬ 
tudinal waves in the column. For small initial disturbances Xq and finite compression 
rates v parameter t*, defined by formula (14), is large. Then using asymptotic formu¬ 
las (13) in (14) yields the explicit dependence of the time required for buckling f* on 
the compression rate: 


f = g(v)(T- 
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where t s = ao/v s . For small initial disturbances Xq and finite velocities v the depen¬ 
dence g(v) is weak. Therefore the time required for buckling is inversely proportional 
to cubic root of the compression rate. 

Consider the contribution of vt* and x* to the critical force §■ For large r* 


formula (14) yields: 


ut*a 0 r*Bi' (r*) 
xl Bi(r*) 


> 1 . 


(16) 


Then the last term in the formula ([9]) can be neglected. 

Substitution of formula (15) into equation (J9| yields the dependence of the critical 
force of the compression rate: 


= l + ^M f—Y 

Fe £e \v s ) 


(17) 


Estimate critical deformation Ee for a column. For Bernoulli-Euler beam with square 
cross-section, thickness h and length l the expression for Ee has the form: 


7T 
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(18) 


Formulas (17) and (18) yield the expression for critical force: 


F* 12g(v)(iyf v y 
F e tt 2 \h) \v s J ’ 


P = 2, 



(19) 


The dependence of g on v and l/h is relatively weak. Therefore critical force and 
compression rate are related by the power law dependence. This result is in a good 
agreement with experimental data BIZ). In paper [7J compression of columns with 
constant rate in a hydraulic machine was investigated. In paper [5J buckling under the 
action a heavy body falling on the end of the column was considered. In both cases the 
dependence of critical force on the compression rate is approximated by formula, similar 
to ( [T7| ). The following values of parameters p, q were obtained: p = 2.0, q = 0.68 [5j 
and p = 1.61, q = 0.62 [7]. 

Estimate the interval of compression rates, where formula ( Jl7| is applicable. The 
derivation is based on the assumption C 2 7^ 0. It is straightforward to show that C 2 = 0 
at some points of the interval Te < —1- Then compression rate should satisfy the 
following condition: 


y_ > 4 

v s l- e e ' 


( 20 ) 


Additionally, the assumption r* > 1 was used, 
yields: 


Using this condition in formula (14) 


v 7r 3 Ai , (rE) 3 

v s > (1 - £e) 4 



Bi(l)3Bi'(l)^. 


( 21 ) 


Thus formula (17) describes the behavior of the system at velocities and initial distur¬ 
bances satisfying conditions (20) and (21). 
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3 Comparison with numerical solution 


The model described above involves significant simplifications. In the present section 
the analytical results are compared with numerical solution of the exact nonlinear 
equation (J4]) . The equation of motion Q is solved numerically using Verlet symplectic 
integration scheme. Initial velocity is equal to zero, initial displacement is equal to xq. 
The following values of dimensionless parameters are used Ct/cl = 1.64-10 -6 , A t/t s = 
10 -2 . The given ratio of stiffnesses corresponds to the column with length/thickness 
ratio about 10 3 (see formula (18)). Velocity belongs to the interval v/v s G [10~ 8 ; 1CT 3 ] 
considered in laboratory experiments [7]. 

In the simulations the longitudinal force is computed using formula (J5]) . Typical 
dependence of the force on deformation for xo/a = 10 -4 and v/v s = 10 _1 °,10 —9 ,10~ 8 
is shown in figure [2} Critical force F* corresponds to the maximum value of the force 
shown in figure [2} Though the velocity is much lower then the “velocity of sound” v s , 



Figure 2: Dependence of longitudinal force on deformation for Xo/a = 1CT 4 , v/v s = 
ICE 8 (dash-dot line), ICC 9 (dashed line), ICC 10 (solid line). Horizontal line corresponds 
to the Euler static force. 


the critical force significantly exceeds Euler static force ([3]). 

The dependence of critical force on the rate of compression v for initial distur¬ 
bances xo/ao = 10~ 3 , ICE 12 is shown in figure [3] It is seen that the dependence is 
accurately approximated by the power law ©■ The power tends to 2/3 as the initial 
disturbance tends to zero. 


4 Conclusions 

Simple one-dimensional model for dynamic buckling of a column under constant rate 
compression is derived. Despite the simplicity, the model qualitatively the influence 
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Figure 3: Dependence of critical buckling force on compression rate for initial dis¬ 
turbances Xo/ao — 10~ 3 (circles), 10~ 12 (diamonds). Points correspond to numerical 
solution, lines correspond to analytical solution (formula (17)). 


of the main parameters (compression rate, length/thicknessration, and initial distur¬ 
bance) on the critical force. In dynamics the critical force buckling exceeds the Euler 
force even at relatively slow compression rates. The reason is that some time is required 
for transition from unstable straight configuration to stable bend configuration. It is 
shown that this time is inversely proportional to cubic root of velocity of column end. 
For small initial disturbances the dependence of critical force on velocity of the end is 
accurately approximated by the power law with exponent approximately equal to 2/3. 
This result is in a good agreement with the results of laboratory experiments 013- 
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